On the optimal form of branching supply and collection networks 
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For the problem of efficiently supplying material to a spatial region from a single source, we 
present a simple scaling argument based on branching network volume minimization that identifies 
limits to the scaling of sink density. We discuss implications for two fundamental and unresolved 
problems in organismal biology and geomorphology: how basal metabolism scales with body size 
for homeotherms and the scaling of drainage basin shape on eroding landscapes. 
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In both natural and man-made systems, branching net- 
works universally facilitate the essential task of supplying 
material from a central source to a widely distributed 
sink population. Branching networks also underlie the 
complementary process of collecting material from many 
sources at a single sink. Such networks typically exhib- 
it structural self-similarity over many orders of magni- 
tude: river networks drain continents [l|, 0, 0|, arterial 
and venal networks move blood between the macroscopic 
heart and microscopic capillaries Q , and trees and plants 
orient leaves in space taking on the roles of both structure 
and transportation. 

We address the following questions regarding supply 
networks: (1) What is the minimum network volume 
required to continually supply material from a source to 
a population of sinks in some spatial region 17?; and (2) 
How does this optimal solution scale if is rescaled allo- 
mctrically? (For convenience, we use the language of dis- 
tribution, i.e., a single source supplying many sinksj Our 
approach is inspired by that of Banavar et al. [y| who 
sought to derive scaling properties of optimal transporta- 
tion networks in isometrically growing regions based on 
a flow rate argument; Banavar et al.'s approach followed 
the seminal work of West et al. who suggested supply 
networks were key to understanding the metabolic limi- 
tations of organisms, and focused on network impedance 
minimization (see @, Q). In contrast to this previ- 
ous work, our treatment is explicitly geometric. We 
also accommodate four other key features: the ambient 
dimension, allometrically growing regions, variable sink 
density, and varying speed of material transportation. 

We consider the problem of network supply for a gen- 
eral class of c?-dimensional spatial regions in a D > c? 
dimensional space. Each region has volume V and 
overall dimensions Li x L2 x ■ ■ ■ x (see Fig. [T^). We 
allow these length scales to scale as Li oc V'^ creating 
families of allometrically similar regions. For isometric 
growth, all dimensions scale uniformly meaning 7i = 1 /d, 
while for allometric growth, we must have at least one of 
the {7^} being different (Fig. [Hd). For the general case 
of allometry, we choose an ordering of {7^} such that the 
length scales are arranged from most dominant to least 




FIG. 1: (a) We consider families of d-dimensional spatial 
regions that scale allometrically with Li oc , and exist in a 
_D- dimensional space where D > d. For the d=D—2 example 
shown, 7max = 71 > 72, and Li grows faster than I/2. We 
require that each spatial region is star-convex, i.e., from at 
least one point all other points are directly observable, and 
the single source must be located at any one of these cen- 
tral points, (b) Distribution (or collection) networks can be 
thought of as a superposition of virtual vessels. In the exam- 
ple shown, the source (circle) supplies material to the three 
sinks (squares) . (c) Allowing virtual vessels to expand as they 
move away from the source captures a potential decrease in 
speed in material fiow. For scaling of branching network form 
to be affected, the radius r of a virtual vessel must scale with 
vessel length s (measured from the sink) as s"*^. 



dominant: 7,„ax = 7i > • • ■ > 7d- 

We assume that isolated sinks are located throughout 
a contiguous spatial region Q (volume V) which contains 
a single source located at 5" = 0. We allow sink den- 
sity to follow p ~ po{V){l + a||a?||)~'» where a is fixed, 
C > 0, and ||.t|| is the distance from the source. When 
the exponent C = 0; P is constant throughout the region 
(as for capillaries in organisms), but remains a function 
of the region's overall volume V. While decreasing sink 
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density > 0) docs not reflect the reality of biolog- 
ical organisms, it is not an unreasonable postulate for 
other supply /collection systems, and may be of use in 
modelling transportation to and from cities. Last, we 
assume each sink draws approximately the same amount 
of material from the source per unit time. The material 
travelling from the source to a specific sink takes up a 
certain volume of the network, and while this volume of 
material may not be coherent away from the sink, we can 
nevertheless imagine separate 'virtual vessels' transport- 
ing material from the source to the sinks (sec Fig. [T})). 
(Only in the smallest, outer branches will virtual vessels 
coincide with physical vessels.) Material flow rate will 
then vary according to changes in the cross-sections of 
these vessels. 

We take the cross-sectional area of these virtual ves- 
sels to be bounded by a fixed upper limit at the sink 
(see Fig.IlJ)). We allow that material speed may increase 
with proximity to the source, meaning these virtual ves- 
sels may taper. If the radius decreases as r-sink(l + cs)~'^, 
where s is the length of the vessel as measured from the 
sink and c is a constant, then the volume of a virtu- 
al vessel grows as f vessel ~ s^~'^'^ for < e < 1/2 and 
I'vessei for e > 1 /2 (we can therefore focus on e = 1/2 
to represent the latter case). We ignore all other possi- 
ble taperings since only an algebraic decay relationship 
between vessel radius and length will affect the scaling 
of overall network volume V^iet- If however, there is a 
minimum virtual vessel radius (i.e., a limit to material 
speed) then vessel volume must grow linearly with length: 

"^vessel ^ ^ ■ 

The overall network volume Vnct is the sum of all virtu- 
al vessel volumes, and is evidently minimized when virtu- 
al vessels travel directly from the source to each sink — the 
extreme case of a star network. While real, large-scale 
distribution networks are branched, many examples are 
close to this limit in terms of path length [13] ■ Minimal 
network volume therefore grows as 

minKotoc/ poil + a\\x\\)-'^\\x\\^-^'dx, (1) 

where < e < 1/2, and we have indicated a spatial region 
scaled to have volume V by ild.oiV). The integral's 
leading order behavior gives the optimal scaling of Vnct 
with V: 

minVnct - PoV'+''"-^'-''-'^\ (2) 

where again po — Pq{V). When e > 1/2, minV^iet ^ 
po^^"''''^"""' ■ The scaling of minimal network volume with 
V is thus governed by sink density p, vessel scaling, and 
the dominant length scale through 7max, and we first 
address the role of the latter. Since for isometric scaling, 
7max = 1/rf, whereas for allometric scaling, 7niax > 1/rf, 
we immediately see that from a scaling perspective, iso- 
mctrically growing regions require less network volume 



than allometrically growing ones, and are in this sense 
more efficiently supplied. Efficiency also increases with 
the dimension d since network volume scales more close- 
ly with overall volume (7max = 1/rf decreases) . Further- 
more, shapes that scale allometrically effectively func- 
tion as lower-dimensional, isometrically scaling objects 
and are therefore less efficiently supplied (the equivalent 
spatial dimension is l/7max)- 

We see from Eq. ([2]) that network volume straightfor- 
wardly increases linearly with pQ. How po in turn scales 
with V depends on the specific system, and in particular 
on whether D ^ d ot D > d. We now specialize our gen- 
eral result for the two cases of blood networks and river 
networks. 

Blood networks [D—d—3]: If material is costly, as 
in the case of blood, then we expect that isometric scaling 
(7max = 1 /d) to be attained by evolution. We take C = 
because capillaries (the sinks) are distributed relatively 
uniformly. Furthermore, since both blood velocity in the 
aorta and in capillaries change little with V during rest- 
ing states (3,[lli, we have e = 0. Lastly, it is well observed 
that the volume of blood scales linearly with organismal 
volume [l^ . Vnet V. (In general, for D — d, wc must 
have Vnet oc as otherwise wc would have the nonsen- 
sical limits of Vnct/V or cx) as oo.) Since 
we already have that Vnet oc poV^'^^^'^, the additional 
constraint Vnct oc V means sink density must decrease 
as volume increases: po (x. V~^l'^ ^ where po for resting 
organisms now refers to the effective or active sink den- 
sity (organisms at rest have a substantial proportion of 
inactive capillaries that are called into use during higher 
rates of overall activity [l3l)- 

It follows that Prest, the average rate of energy use 
in a resting state (basal power), which is proportional 
to the number of active sinks in fi, can at best scale 

we have 

assumed that V scales as body mass M . For three dimen- 
sional organisms, we therefore have 

P,est K Af2/3. (3) 

If organism shapes obey instead an allometric scaling 
then power scales more slowly as Prest cx Af with 
7max > 1/3, contrary to the McMahon's theory of elastic 
similarity |14l |. 

We note that in detail, blood networks do not appear 
to have universal forms [l5| . showing substantial varia- 
tion in branching structure across and within species; we 
therefore argue that it is only the system level that mat- 
ters and that branching networks need only approximate 
star networks. 

Crucially, the scaling law of Eq. ([3]) balances with the 
standard one based on organismal surface area S. For 
homeothermic organisms, who must constantly balance 
heat loss to maintain a steady internal temperature, we 
have that Prest oc S (due primarily to radiation but also 
convection [lij). For isometrically scaling organisms. 
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S oc Af^/^, and this is well supported empirically [l^ . 
Moreover, it is easy to show that only isomctrically grow- 
ing shapes balance Prcst since the d — 1 dimensional sur- 
face area of a growing region scales as T/i-'^-"-, and 
therefore 7max = 7min = Thus, the most efficient 

network in terms of minimal volume is also the one that 
precisely balances radiative heat loss. 

Our seemingly reasonable result, which was empirically 
observed over a century ago by Rubner , runs counter 
to nearl y e ighty years of reports that Prest oc M^/"'. 
Kleiber [l8| first suggested the exponent might be 3/4 
in the 1930's after measuring a value of 0.76 for 13 mam- 
mals (his practical reason for choosing 3/4 was to sim- 
plify slide rule calculations In the decades follow- 
ing, a general but not universal consensus on a "3/4- 
law of metabolism" was reached 0, [20l. j2ll The issue 
remains controversial both theoretically [sl. |22|. [23l . [23 |. 
and empirically: some recent statistical analyses have 
shown that a 2/3 exponent is well supported by large data 
sets for warm-blooded organisms (both birds and mam- 
mals) [1, [2^ [2^, [l^l while others have found evidence 
in favor of a 3/4 exponent or no simple scaling relation- 
ship d, M, S M, l3l| : and for cold-blooded organisms, 
plants, and invertebrates, a much broader range of expo- 
nent values have been measured empirically and predict- 
ed from theory 

Perhaps the most important aspect of the 2/3 vs. 3/4 
debate is that from an optimization point of view, the 
lower the scaling of resting metabolism the better. A 
scaling of M^^^ (or any power exceeding 2/3) would 
point to either a fundamental scaling limitation for warm- 
blooded organisms, or to the existence of a cost oth- 
er than volume minimization, such as impedance 
Where a higher exponent would be desirable is in the 
scaling of maximal power Pmax, which is unsustainable 
and depends on stored en ergy , and indeed, Pmax scales 
almost linearly with mass [33[ . 

River networks [D—3 > d—2]: The patterns of 
large-scale river networks have long drawn scientific 
interest, naturally from hydrologists and geomorpholo- 
gists [l|, but also from statistical physicists seeking 
evidence of universality in nature [s^, [s^, [H, [s^ . In our 
framework, river networks are collection systems: water 
flows from many sources (channel heads [40|) to a sin- 
gle sink, the outlet of the network's main stream. The 
description of river network geometry has often focused 
on Hack's law [4lj which relates the area a of a drainage 
basin to the length I of its longest stream: I oc a'\ Var- 
ious studies of small-scale basins, starting with Hack's 
initial work [1, [4l[, have suggested that the 'Hack expo- 
nent' h exceeds 1/2, indicating an anomalous allometric 
scaling of basin shape wherein large basins are relatively 
long and thin compared to smaller ones. Indeed, for suf- 
ficiently small, homogeneous landscapes, there may be a 
universal value oi h > 1/2, yet to be fully understood 
theoretically [sl, 37, 3^. However, the most comprehen- 
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FIG. 2: The scaling of network volume versus basin area for 
four continental-scale river networks The solid line indicates 
a scaling of 3/2. Network volume is estimated for an idealized 
steady-state, uniform rainfall condition normalized such that 
V^iGt = JZij > where aij is the area of the basin draining 
into the ij'th cell on a coarse-grained version of a landscape, 
and the data is binned in log space [^. 



sive data sets robustly show that ft, = 1/2 for large-scale 
networks; in particular, Montgomery and Dietrich (40| 
found that over twelve orders of magnitude variation in 
basin area, / oc a^'"*^ (their data set mixed both Euclidean 
overall basin length L and main stream length Z; gener- 
ally, / oc with /3 very close to, if not equal to, unity). 

The empirical observation that h = \/2 accords with 
our result that with respect to network volume mini- 
mization, isomctrically growing regions arc most efficient; 
with our optimality argument, this becomes a stronger 
statement than appealing only to dimensional analysis. 

Beyond isometry, we have the scaling of network vol- 
ume to consider. We now take po to be constant and 
again set Q = Q, meaning we assume that, when aver- 
aged over time, rain falls approximately uniformly across 
a landscape. We also assume e = for the continental 
scale networks we examine below. In contrast to the case 
of cardiovascular networks, the constraint that network 
volume must scale as overall volume (or basin area a) 
cannot apply to river networks, since Eq. ^ now gives 
Kot oc piiV^^^/'^ oc y^/^ = a^/^. The reason is simple: 
river networks lie on a (i = 2 dimensional surface embed- 
ded in D = 3 dimensions, and the presence of a third 
dimension allows the total water in the network to grow 
faster than if the embedding dimension was D ^ 2. 

In Fig. [2] we show that Vnet scales as a^^^ for four 
continental-scale networks: the Mississippi, the Amazon, 
the Congo, and the Nile [i^. The scaling is robust and 
holds over 6 and 10 orders of magnitude in area and 
network volume respectively, indicating that self-similar 
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drainage basins most cfRcicntly drain large-scale land- 
scapes. 

While we have argued that the optimal scaling corre- 
sponds to the isometric case of /i = 1/2, we can directly 
see this connection arising from known scaling laws of 
river networks 3^, 3^. As put forward in (37| . network 
volume scales with basin area as Vnct oc a^~^^, showing 
that for the optimal case, we indeed have h = 1/2. Note 
that only by considering allometrically growing regions, 
as we do here, can this connection be formally made. 

We acknowledge that a stronger optimization may be 
at work for eroding landscapes. In particular, previ- 
ous theoretical analyses suggest that in terms of ener- 
gy minimization, landscapes may reach local, dynami- 
cally accessible minima 37, [s^, though sufficient self- 
averaging inherent in larger landscapes may render such 
minima indistinguishable from the global one. Neverthe- 
less, our argument and data analysis show that for large- 
scale networks on heterogeneous landscapes, far exceed- 
ing the typical correlation length for precipitation events, 
a kind of optimal volume minimization is achieved. 

With suitable modifications, our findings may be found 
relevant to other systems, in particular to plants (when 
seen as two connected branching networks) , as well as to 
the scaling limits of episodic movement such as the trans- 
portation of people in and out of city centers. (Assuming 
constant flow as we do here, then a decay of sink density 
with C = 1 follows for growing cities where transporta- 
tion remains unchanged with city size.) While we have 
demonstrated that empirical evidence supports a geomet- 
ric optimality for two kinds of large-scale natural branch- 
ing networks, not all systems will be optimal or may be 
optimized. For example, if allometric scaling of organis- 
mal shape is demanded by some other constraints (e.g., 
due to the effects of gravity) , then blood volume will be 
forced to obey a poorer scaling with overall volume Q . 
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